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There have been very few measurements of the total cross section for np scattering below 500 
keV. In order to differentiate among NN potential models, improved cross section data between 20 
and 600 keV are required. We measured the np and nC total cross sections in this energy region by 
transmission; a collimated neutron beam was passed through CH2 and C samples and transmitted 
neutrons were detected by a BC-501A liquid scintillator. Cross sections were obtained with a 
precision of 1.1-2.0% between 150 and 800 keV using ratios of normalized neutron yields measured 
with and without the scattering samples in the beam. In energy regions where they overlap, the 
present results are consistent with existing precision measurements, and fill in a significant gap in 
the data between E n — 150 and 500 keV. 



I. INTRODUCTION 

Nucleon-nucleon interactions and NN po- 
tential models are an important representation 
of the strong interaction and a component of 
the theory of nuclei. Measurements of the two- 
body interaction in pp and np reactions have 
been applied to studies of the structure and 
dynamics of light nuclei as well as testing the 
charge, spin, and isospin dependence of the 
strong nuclear force [l| . 

Additionally, np scattering data are appli- 
cable to the detection of neutrons in organic 
scintillators below 2 MeV, where the primary 
mechanism for depositing energy in the scintil- 
lator is np elastic scattering. Therefore, deter- 
mining the efficiency of a neutron detector us- 
ing either Monte Carlo or analytical methods 
requires accurate knowledge of the np elastic 
cross section. Since the np — > dry cross sec- 
tion is typically 5 orders of magnitude smaller 
than the clastic cross section in this energy 
range Q, measurement of the total np cross 
section effectively yields the elastic cross sec- 
tion. 
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However, the relevant data iH plotted 
with the ENDf/B-VII.l tabulation j|] in fig- 
ure [TJ reveal a lack of precise measurements 
between 150 and 500 keV. Previous measure- 
ments all employed the method of neutron 
transmission to determine the total cross sec- 
tion, using polyethylene and carbon samples to 
account for the carbon contribution and obtain 
the np cross section. We focused on the range 
between 150 and 800 keV, covering the region 
below 500 keV where there are very limited 
data. 



II. EXPERIMENTAL PROCEDURE 

We measured the total cross section for neu- 
tron scattering from polyethylene and pure 
carbon samples using the method of neutron 
transmission. The neutron beam were pro- 
duced via the 7 Li(p, n) 7 Be reaction, with a 
1.875 MHz pulsed proton beam produced by 
the Van de Graaff accelerator located at the 
University of Kentucky. The Q-value for this 
reaction is -1.644 McV [?|]. The neutron pro- 
duction target was composed of 20 kA thick 
LiF on a tantalum backing; this thickness was 
chosen such that 2.25 MeV protons would lose 
approximately 50 keV before exiting the tar- 
get. Using a thin target limited the spread 
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FIG. 1. Existing np total cross section measurements for neutron energies from 150 to 800 keV, measured 
by (squares) Frisch (circles) Bailey et. al. [J, (triangles) Bretscher et. al. [f|, and( diamonds) Clement 
et. al. 0| plotted with the ENDf/B-VII.l tabulation (solid line) from Hale et. al. @|. 
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FIG. 2. Detector and shielding configuration. 

of neutron energies and reduced the back- 
grounds. 

The neutrons were collimated by a copper 
shield, with a 6.35 cm opening and 50.8 cm 
thickness. The beam was additionally defined 
by a wax collimator with an opening matched 
to that of the copper shield and tapered down 
to a 2.2 cm opening. The sample was placed in 
the beamline past this collimation. The neu- 
trons were detected by a 13.7 cm diameter BC- 
501A liquid scintillator, which was housed in a 
lead and wax shield, with a 4.4 cm lead lining 
in the internal aperture and a 7 cm wax lin- 
ing around the lead. The experimental setup 
is shown in figure [5] The geometry of the col- 
limation resulted in a beam size smaller than 
both the samples and the neutron detector. 
This was tested by performing several runs 



where a sample was shifted by several cen- 
timeters off its centered position in the beam. 
These shifts did not produce any changes in 
the neutron yields, indicating that the neutron 
beam is confined to the central region of the 
samples. Previous experiments @,|4|], in which 
the beam profile was larger than the detector, 
required an additional correction due to neu- 
trons produced at angles outside the detector 
acceptance which scattered into the detector. 
In our geometry, this correction was not re- 
quired. 

The bombarding proton energies ranged 
from 2.00 to 2.55 McV in 50 keV increments, 
yielding neutrons from 150 to 800 keV. This 
spacing produced minimal overlap of neutron 
energies, especially above a proton energy of 
2.25 MeV. This resulted in an increase in the 
statistical uncertainty near neutron energies 
of 650, 700, and 750 keV. For proton ener- 
gies above 2.25 McV (neutron energies above 
450 keV), a total of four hours of data were 
taken at each proton energy. For proton ener- 
gies below 2.25 MeV (neutron energies below 
450 keV) a total of eight hours of data were 
taken at each proton energy. These additional 
statistics were collected to compensate for the 
decreasing 7 Li(p, n) 7 Bc cross section and the 
increasing np total cross section. 

Four polyethylene and three carbon samples 
were used, in addition to the empty position 
used to determine the yield with no sample in 
the beam. A caliper was used to measure each 



Material 


Nominal Thickness (cm) 


t (g/cm ) 


CH 2 


0.5 


0.479 ± 0.003 


CH 2 


1.0 


0.959 ± 0.006 


CH 2 


2.0 


1.93 ±0.01 


CH 2 


3.0 


2.91 ±0.02 


C 


3.1 


4.89 ± 0.02 


c 


4.6 


7.24 ± 0.03 


c 


6.1 


9.42 ± 0.04 



TABLE I. CH 2 and C sample thicknesses. 



target dimension with an uncertainty of 0.25 
mm. The mass was determined with a scale 
to an accuracy of 10 mg. The polyethylene 
samples were rectangular prisms, with nomi- 
nal length and width of 7.6 and 5.1 cm, with 
a thickness varying from 0.5 to 3.0 cm. The 
carbon targets were cylindrical, with a diame- 
ter of 4.8 cm and a thickness varying from 3.1 
to 6.1 cm. The nominal thickness and areal 
density of each sample are given in table |U 
Approximately 30% of the transverse area of 
each sample was exposed to the neutron beam. 

The samples were mounted on a remotely- 
controlled wheel which could position up to 
twelve samples in the neutron beam. Dur- 
ing each two hour run, the wheel would au- 
tomatically switch which target was in the 
beam at pre-programmed intervals, resulting 
in a series of irradiations between one and 
seven minutes long. Feedback signals from 
the wheel indicated when samples were being 
moved, and which sample was in the beam 
at a given moment. By using these short 
intervals, the yields from each irradiation of 
the samples could be compared under similar 
beam conditions, eliminating potential varia- 
tions due to beam current fluctuations or the 
condition of the LiF target. At each of the 
twelve positions, the samples were mounted 
onto a smaller wheel. When in the beam posi- 
tion, a second motor rotated this smaller wheel 
so as to average over possible non-uniformities 
in thickness. 

Additionally, a sulfur sample was used to 
calibrate the neutron energy. With its multi- 
ple resonances in our energy range, we can use 
the observed decreases in yield and the max- 
imum neutron energy for each proton beam 
energy to calibrate the absolute time of flight. 
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FIG. 3. Neutron energy spectrum with a sulfur 
sample for E p = 2.00 MeV. This corresponds to 
a maximum neutron energy of 230 keV, and there 
is a visible decrease in counts associated with the 
sulfur resonance at E n — 203 keV. 



The neutron energy spectrum with a sulfur 
sample for E p = 2.00 MeV is shown in fig- 
ure [21 

The data were recorded event by event, us- 
ing NIM and CAM AC electronics. The dead 
time ranged from 5-40%, depending on the 
event rate. 



III. ANALYSIS 

We used pulse shape discrimination in the 
neutron detector to separate neutrons and 
7-rays, thus significantly reducing the back- 
ground. The pulses from the liquid scintillator 
were split and were recorded in separate ADCs 
using a long gate (500 ns) and a short gate (100 
ns). The short-gated pulse height versus the 
long-gated pulse height is shown in figure 2] 
The upper and lower broad bands correspond 
to 7-rays and neutrons, respectively. They 
merge at small pulse heights, where the differ- 
ence in the long and short portions of the pulse 
is too small to separate, but using a conserva- 
tive cut allows the elimination of all higher en- 
ergy 7-rays. Figure [5] shows the neutron time 
of flight spectrum for E n = 450 to 500 keV 
with no conditions (solid) and the pulse shape 
discrimination condition (dashed). The back- 
ground outside the neutron peak (270-300 ns) 
is reduced by 75% while the neutron yield is 
unaffected. Neutron yields were determined 
in 10 keV bins from 150 to 800 keV, and the 
integrated live current was recorded for each 
sample irradiation in order to normalize the 
neutron yields. 

To extract a cross section from the yields 
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FIG. 4. (Color online) Neutron detector ADC out- 
puts: short-gated pulse height versus long-gated 
pulse height. The upper broad band is 7-rays, the 
lower band is neutrons. 
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FIG. 5. Neutron time of flight for E n = 450 to 500 
keV with no conditions (solid line) and the pulse 
shape discrimination condition (dashed line) . The 
7-fiash is visible at 10 ns, with the neutron peak 
visible at 270-300 ns. With the pulse shape dis- 
crimination condition, the background is reduced 
by 75% outside the neutron peak. 



from each sample, we normalized the neutron 
yields by the proton beam current, 
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where t p and r c are the sample thicknesses 
from table HI Y is the neutron yield, and Q is 
the integrated live-time current. The carbon 
contribution of the polyethylene targets is sub- 
tracted, and the uncertainty in the nC cross 
section is included in the total uncertainty in 
the np cross section. Due to the long lengths 
of the CH2 chains, the variation of the hy- 
drogen to carbon ratio from 2 is insignificant 
compared to the other systematic uncertain- 
ties. These results do not require knowledge 



of the neutron detector efficiency nor the ab- 
solute beam flux. The proton beam current 
was integrated for each irradiation, and the 
dead time of the data acquisition system was 
precisely measured as a function of the event 
rate in order to determine the total live-time 
current. 

The backgrounds were divided into two cat- 
egories: the room background, which is inde- 
pendent of the time of flight, and the sample 
background, which is due to neutrons which 
scatter from the sample but are still detected. 
The room background can be measured from 
the data by fitting to the constant background 
outside the neutron time-of-flight peak. The 
sample background was determined by extrap- 
olating to zero thickness. For each target, the 
cross sections measured at each energy were 
combined using a weighted average, in order 
to produce a single average cross section for 
the entire energy range. These average cross 
sections were then plotted versus target thick- 
ness and fitted with a linear function. The 
intercept of this linear function gave the aver- 
age cross section at zero target thickness. By 
taking the ratio of the average cross section at 
zero thickness and the average cross section for 
each sample, we determine the correction for 
the background due to target thickness. These 
corrections ranged from 0.3% for the thinnest 
sample to 1.3% for the thickest sample. The 
corrected cross sections for each energy bin 
were then combined using a weighted average. 



IV. RESULTS 

Our results are shown in figure [5^ for the 
np total cross section and in figure Uh- for the 
nC total cross section. Existing np data 
and the ENDf/B-VII.l tabulation are also 
shown in figure^, and existing nC data [il-tioj 
and the ENDf/HE-VI [U tabulation are also 
shown in figure [7^,. Figures |6Jd and [TJd show 
the difference between our experimental cross 
sections and the tabulations. Error bars are 
suppressed in the previous nC measurements 
for clarity in the plot; the uncertainty in these 
measurements is between 2.7% and 5%. 

The total uncertainty in each of our mea- 
surements is in the range 1.1-2.0%. The sys- 
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FIG. 6. (a) Results for (filled circles) the np total cross section measurements for neutron energies from 
150 to 800 keV, plotted with ENDf/B-VII.l tabulation (solid line) from Hale et. al. @] and previously 
measured np cross sections from 150 to 500 keV, measured by (squares) Frisch (open circles) Bailey 
et. al. [U, (triangles) Bretscher et. al. and (diamonds) Clement et. al. [f|. (b) The difference between 
our experimental cross section and the ENDf/B-VII.l tabulation. Error bars include the statistical and 
systematic errors added in quadrature, along with the contribution due to subtracting the experimental 
nC cross section. 



tematic uncertainties included the measured 
masses and dimensions of the samples, listed 
in table HI contributing 0.35%; the characteri- 
zation of the beam and system deadtime, con- 
tributing 0.5%; and the background determi- 
nation, contributing 0.5%. The statistical un- 
certainty was of the order of 0.4%. The uncer- 
tainties due to the target dimensions applied 
uniformly to all data points, while the statisti- 
cal, beam, system dead-time, and background 
uncertainties were determined for each sample 
irradiation. 



V. DISCUSSION 

In order to connect our measurements to 
NN potential models, we parameterize the s- 
wave neutron-proton clastic scattering cross 
section, 

3 1 

° = 4°* + I' 7 "' ^ ' 

in terms of the triplet and singlet scatter- 
ing lengths a t s and e nerg y-dependent effective 
ranges p M (0,T): [H,[lf 



(a-; 1 - y d (0,T)p*Y+p*y W 

where the subscript d represents either t or 
s, pd(0,T) is the energy-dependent effective 
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FIG. 7. (filled circles) Results for the nC total cross section measurements for neutron energies from 
150 to 800 keV, plotted with the ENDf/HE-VI tabulation [rfl | and existing data, measured by (squares) 
Huddleston et. al. (open circles) Wilenzick et. al. and (triangles) Uttley et. al. Error bars 
are suppressed in the previous nC measurements for clarity in the plot, (b) The difference between 
our experimental cross section and the ENDf/HE-VI tabulation. Error bars include the statistical and 
systematic errors added in quadrature. 



range, and T and p are the center of mass 
kinetic energ y a nd momentum, respectively. 
Hackenburg [13( has recently revisited the 
problem of determining the zero-energy cross 
section and effective-range theory (ERT) pa- 
rameters from data, including a considera- 
tion of the correlation between the singlet and 
triplet effective range. Measurements of the 
zero-energy np cross section, crn [HI, EH , and 
the parahydrogen coherent scattering length, 
a c [1103, given by 

(J = 7r(3fl2 + a 2), ^ 

a c = § a t + 5 a si 

were included in the fit due to correlation be- 
tween a s and a t - The ERT parameters result- 
ing from the fit were p t (0, 0) = 1.718 ± 0.025 



fm and p s (0,0) = 2.696 ± 0.059 fm. Also, 
the zero-energy shape dependence parameter 
Ar t = -0.025 ± 0.025 fm. Hackenburg's anal- 
ysis also concluded that additional total cross 
section measurements for np scattering be- 
tween 20 and 600 keV were required in order 
to use the parameters of effective range theory 
to differentiate among NN potential models. 

It was hoped that the present results in the 
150-600 keV energy range would reduce the 
uncertainties in the parameters determined 
from the fit to the ERT expression (Eq. 0]) 
for the cross section. Following Hackenburg's 
method, we included our data in the fit. 
The resulting parameters were consistent with 
those determined by Hackenburg. The uncer- 
tainties were not significantly reduced, as our 
measurements were not precise enough to fur- 
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ther constrain the parameters. According to 
Hackenberg, a 0.004% precision measurement 
at 130 keV is required to reduce the uncer- 
tainty in #(0,0) and Ar t to 0.001 fm [3. By 
instead measuring the energy dependence of 
the cross section, we estimate that a precision 
of 0.5% across the 150 to 800 keV range cur- 
rently measured is required to decrease these 
uncertainties, and 0.1% measurements across 
this range will decrease the uncertainty to 0.01 
fm. 



VI. CONCLUSION 

Our measurement of the np total scattering 
cross section has filled in a large gap in the to- 
tal cross section measurement below neutron 
energies of 500 keV. By measuring ratios of 
transmitted events with and without the sam- 
ples in the beam, we were able to determine 
the cross section independently of the neutron 
detection efficiency. Both the np and nC to- 
tal cross sections are consistent with previous 
measurements, and the nC results show signif- 
icantly decreased scatter and uncertainty. 

There is a slight systematic discrepancy be- 
tween our measurements and the ENDf tab- 
ulations, with a —18 mb shift downward in 
the nC cross section, and a +136 mb shift up- 
ward in the np cross section. This discrepancy 
would not have been visible in the previous nC 



data. The previous np data in this range also 
show a slight increase relative to the theoreti- 
cal curve. 

However, when the present data are used to 
try to improve the fit of effective range theory, 
we find that our measurements are not precise 
enough to increase the precision on the result- 
ing parameters. The effective range parame- 
ters pd are most sensitive to measurements in 
the 20-600 keV region, but more precise data 
at higher energies have already constrained the 
parameters more narrowly than the precision 
of current measurements in this region. The 
discrepancy noted previously is comparable to 
the current uncertainty, and so does not signif- 
icantly impact the determination of the ERT 
parameters. To be comparable to the current 
uncertainty in pt(0, 0) and Ar t , 0.5% preci- 
sion would be required across the energy range 
measured here, 150 and 800 keV. This preci- 
sion is attainable with increased statistics and 
better characterization of the beam, both of 
which are possible with the current configura- 
tion. 
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